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Abstract 

For very large values of tan j3, the charged Higgs boson pair production via bb annihi- 
lation can proceed dominantly at the Large Hadron Collider (LHC) . We calculated the 
cross sections of the charged Higgs boson pair production via subprocess bb — > H + H~ 
at the LHC including the next-to-leading order (NLO) QCD corrections in the minimal 
supersymmetric standard model (MSSM). We find that the NLO QCD corrections can 
significantly reduce the dependence of the cross sections on the renormalization and fac- 
torization scales. 
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I. Introduction 



One of the most important missions of future high-energy experiments is to search for scalar 
Higgs bosons and explore the electroweak symmetry breaking mechanism. In the standard 
model (SM)jJ, one doublet of complex scalar fields is needed to spontaneously break the 
symmetry, leading to a single neutral Higgs boson h°. The minimal supersymmetric standard 
model (MSSM) [2] is one of the most attractive extensions of the SM. The MSSM requires the 
existence of two doublets of Higgs fields to cancel anomalies and to give masses separately to 
up and down-type fermions. The MSSM predicts two CP-even neutral Higgs bosons h°,H°, 
a pseudoscalar A Higgs boson and a pair of charged scalar particles H^. At the tree level, 
the MSSM Higgs sector has two free parameters: tan/3 = vz/vi, the ratio of the vacuum 
expectation values of the two Higgs doublets and a Higgs boson mass which is taken to be 
mjj± in this paper. 

The discovery of the would be a clear signal for the existence of physics beyond the 
SM with a strong hint towards supersymmetry. The CERN large hadron collider (LHC), with 
yfs = 14 TeV and a luminosity of 100 fb^ 1 per year, will be a wonderful tool for looking for 
new physics. At the LHC, the light charged Higgs boson can be produced from the top quark 
decays t —* b + -ff + [3]. Heavy charged Higgs boson is mainly produced via the processes 
gb — > gg — * tbH + ^\ and qb — * q'bH + ®. Moreover, single charged Higgs boson 

production associated with a W boson, via tree-level bb annihilation and one-loop gg fusion, 
has been proposed and analyzed in Ref.[7]. 

At the LHC, the charged Higgs boson also can be produced in pair production mode. There 
are three important H + H~ production channels: (i) qq — > H + H~, where q = u,d,c,s,b. 
(via Drell-Yan process, where a photon and a Z-boson are exchanged in the s-channel. In 
the case of q = 6, there are additional Feynman diagrams involving h° and H° in the s- 
channel and the top quark in the t-channel) ,8_ . For very large values of tan /3, due to 
the large contributions from the additional diagrams, the H + H~ production can proceed 
dominantly via bb annihilation 9 . (ii) gg — * H + H~ (via quarks and squarks loop) • 
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(iii) qq — ► qqH + H~ (via vector boson fusion) jllj . 

In the subprocess bb — ► H + H~ , the initial state bottom quarks arise from a gluon in the 
proton splitting into a collinear 56-pair, parameterized in terms of bottom quark distribution 
functions. On the other hand, the 'twin' process gg — ► bbH + H~ using gluon density has 
been studied at LO in Ref.|12j. It is pointed out that the use of the b-quark density may 
overestimate the inclusive cross section due to crude approximations in the kinematics |13| . 
However, it is suggested that the bottom quark parton approximation maybe valid by choosing 
appropriate factorization scale ^1] . Following the suggestions in Ref . |13l I14j , we can analyze 
the transverse momentum distribution of the b-quarks in the process gg — > bbH + H~ as 
shown in Fig. 2 of Ref.|12j. The most suitable factorization scale for bb — > H + H~ is of order 
mji± /5 ~ mu±/A which is much smaller than the usually used scale m#±. 

In this paper, we study the process pp —> bb — > H + H~ + X at the LHC with very large 
tan/3 in the MSSM. The NLO QCD corrections are calculated. The paper is organized as 
follow: In section 2, we discuss the LO results of the subprocess bb — ► H + H~ . In section 3, 
we present the calculations of the NLO QCD corrections. In section 4, the numerical results, 
discussions and conclusions are presented. 

II. The Leading Order Cross Section 

The Feynman diagrams for the subprocess b{p\)b(jp2) — * H + {kz)H~ (k<±) at the LO are shown 
in Fig.l, where p\ 2 an d £3,4 represent the four- momenta of the incoming partons and the 
outgoing particles respectively. 




b 

WO M) (e) 

Figure 1: The tree level Feynman diagrams for bb — > H + H~ . 
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We divide the tree-level amplitude into two parts, 

M° = M { Q s) +M^\ (2.1) 

where and represent the amplitudes arising from the s-channel diagrams shown 

in Fig.l(a,b,c,d) and the t-channel diagram shown in Fig. 1(e) respectively. The explicit 
expressions for the amplitudes and can be written as 

t — m 2 

a^( s ) \f 9H+H-h ( >9h%b , 9h+h-h^9hHI Qbe 2 (h ~ h) 

M — IV{P2)[ 5 1 ' 2 ' 

s — m^ s — m z H0 s 

~ 2c 2 s 2 J ( } 

where s = (pi + P2) 2 , t = (pi — k^) 2 and u = (p\ — k^) 2 are the usual Mandelstam variables. 

Pl,r = (It 75)/2j = sin6* w , c w = cos6 w ,Qb = —1/3. The couplings are defined below, 

(1) iem t ( 2 ) _ ierrib(n r ) tan/3 

~ ^ mw /s w tan/3' ~ y/2m w s w ' 

iemb(fi r ) sin a iemb(p r ) cos a 

= 2m w s u , cos/3' = ~ 2m wSa) cos/3' 

zemvK , . ,„ . cos(2/3) sin (a + /3) , 
9H+H-h° = (sm(/3-a) + — 2 ), 

S w 

iem w cos(2/3) cos(a + /3) 
9h+h-h° = (cos(/3-a) ), (2.3) 

S w 

where a is the mixing angle which leads to the physical Higgs eigenstates h° and H°. mb{/J, r ) 
is MS running mass of the bottom quark. We neglected the bottom quark mass during our 
calculation except in the Yukawa couplings. 

Then the LO cross section for the subprocess bb — > H + H" is obtained by using the 
following formula: 

a o {ijb l^ H+H - ) = __ dt^|M°| 2 , (2.4) 

where i max ,min = ( m H ± ~ 5*) =t \ \js 2 — 4m^-±s. The summation is taken over the spins and 
colors of initial and final states, and the bar over the summation denotes averaging over the 
spins and colors of initial partons. 



III. NLO QCD Corrections 

The NLO QCD corrections to pp -> bb -> H + H~ + V in the MSSM can be separated into 
two parts: the virtual corrections arising from one loop diagrams and the real corrections. 

III..1 Virtual One-loop Corrections 

The virtual one-loop diagrams of the subprocess bb — > H + H~ in the MSSM consist of self- 
energy, vertex and box diagrams which are depicted in Figs. 2-3. Fig.2 shows the diagrams 
of the SM-like QCD corrections arising from quark and gluon loops, and Fig.3 shows the 
diagrams of the so called 'pure' SUSY QCD corrections arising from squark and gluino loops. 
There exist both ultraviolet (UV) and soft/collinear infrared(IR) singularities in the amplitude 
for the SM-like diagrams shown in Fig.2. The amplitude for the 'pure' SUSY QCD diagrams 
(Fig.3) contains only UV singularities. In our calculation, we adopt the 't Hooft-Feynman 
gauge and all the divergences are regularized by using dimensional regularization method in 
d = 4 — 2e dimensions. 



b H+ b H+ 




(g) (h) (i) 

Figure 2: The one-loop Feynman diagrams of the SM-like QCD corrections for bb — > H + H~ 
subprocess. 

In order to remove the UV divergences, we need to renormalize the wave functions of 
the external fields and the Yukawa couplings of h° — b — b, H° — b — b and H~ —t — b. We 
renormalize the top quark mass in the on-mass-shell (OS) scheme. For the renormalization of 
the bottom quark mass in the Yukawa couplings, we employ the modified minimal subtraction 
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Figure 3: The one-loop Feynman diagrams of the 'pure' SUSY QCD corrections for bb 
H + H~ subprocess. 

(MS) scheme. The relevant renormalization constants in this work can be expressed as 
6m t 



m t 



^ m t \(QCD) + ^ m t ySQCD) 
mt mt 



a. 



--^C F {(AB + 2B 1 )(mt,mt,0) - 1} 



a s 
Air 



C F {J2[Bi 



m~„ 



i=i 



m t 



9 sin(2^)(-irB ](m?,m|,m?)}, 



5Z b L 



6Z« 



5rrib 



= (5Zt) {QCD) + (5Zt) {SQCD) 



a s 
4tt 



C F B (0,0,0) + ^-C F [Bi(0, m?,m? ) cos 2 6» 5 + #i(0, m?, m?J sin 2 6» 5 



bi' 



4n 



2ir 



-^-C F B (0,0,0) + -^-C F Bi(0, m|, ) sin 2 0g + Bi(0,m?,m?J cos 2 ^ 



m 6 



47T 47T 



(3-1) 



where A = ^ — 7^ + ln(47r) and C F = |. In above equations we divide the renormalization 
constants into two parts, one arises from the one-loop diagrams involving quark and gluon, 
the other comes from the loops involving squark and gluino. 

The virtual corrections to the cross section for the subprocess bb — > H + H~ can be written 

as 



a v (s, bb — > H + H~) 



1 



167TS 2 



dt 2ReY}{M V ) ] M\ 



(3.2) 



where t maX] min = { m H ± ~ y ^ 2 ~~ 4m^ ± s, and the summation with bar over head means 
the same operation as that appeared in Ea. (|2.4|) . M v is the renormalized amplitude for virtual 
one- loop corrections. After renormalization procedure, a v is UV-finite. Nevertheless, it still 
contains the soft /collinear IR singularities 



where 



A\ = -2C F , A\ = -3C F . (3.4) 

The soft divergences will be cancelled by adding the real gluon emission corrections. The 
remaining collinear divergences can be absorbed into the parton distribution functions, which 
will be discussed in the following subsections. 

III. .2 Real Gluon Emission Corrections of bb — > H + H~ + g 

The real gluon emission subprocess bb — > H + H~+g (shown in Fig. 4) presents 0(a s ) correction 
to bb — > H + H~. It also gives the IR singularities which cancel the analogous singularities 
arising from the one-loop level virtual corrections mentioned in the above subsection. These 
singularities can be either of soft or collinear nature and can be conveniently isolated by slicing 
the phase space of subprocess bb — > H + H~+g into different regions defined by suitable cutoffs, 
a method which goes under the general name of the phase space slicing method fPPS)[T5]. 
We denote this 2^3 subprocess as 

b(p 1 ) + b(p 2 )^H + (k 3 ) + H-(k 4 )+g(k 5 ) (3.5) 

and calculate the cross section by using the two cutoff phase space slicing method [THj. We 
define the Lorentz invariants 

s = (p x + p 2 ) 2 , t = (pi - k 3 ) 2 , u = (pi - /c 4 ) 2 , 
ii5 = (Pi ~ h) 2 , i 2 5 = (j>2 ~ h) 2 (3.6) 
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Figure 4: The Feynman diagrams for bb — > H + H + g. 

and describe this method briefly as follows: Firstly, by introducing an arbitrary small soft 
cutoff 5 S we separate the 2^3 phase space into two regions, according to whether the energy 
of the emitted gluon is soft, i.e. E§ < 5 s V§/2, or hard, i.e. E5 > S s V§/2. The partonic real 
cross section can be written as 

o*(bb -> H + H~g) = o s g (bb -» H + H~ g) + of (bb -► H + H~g), (3.7) 

where Og is obtained by integrating over the soft region of the emitted gluon phase space. Og 
contains all the soft IR singularities. Secondly, to isolate the remaining collinear singularities 
from of, we further decompose of into a sum of hard collinear (HC) and hard non-collinear 
(HC) terms by introducing another cutoff 5 C named collinear cutoff 

of {bb -> H+H-g) = &f G (bb -> H + H~g) + af^ibb -» H + H~g). (3.8) 

The HC regions of the phase space are those where any one of the Lorentz invariants t\$ , £25 
becomes smaller in magnitude than 6 c s, while at the same time the emitted gluon remains 
hard. <r^ c contains the collinear divergences. In the soft and HC region, a~j and <r^ c can be 
obtained by performing the phase space integration in d-dimension analytically. In the HC 
region, <7^ c is finite and can be evaluated in four dimensions using standard Monte Carlo 
techniques 17 1- The cross sections, dg , o"^ c and o"^ 10 , depend on the two arbitrary parameters, 
5 S and S c . However, in the total real gluon emission hadronic cross section of after mass 
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factorization, the dependence on these arbitrary cutoffs cancels, as will be explicitly shown 
in Sec. 4. This constitutes an important check of our calculation. 
The differential cross section in the soft region is given as 

~a s T(l - e) / 



with 



da a g = da v 



2vrr(l -2e) 



(3. 



A$ = 2C F , Af = -AC F ln8, 



si 



The differential cross section da^ c can be written as 



Ag = AC F ln 2 5 s . 



(3.10) 



da 



HC 



da 



s.0 



a s T(l - e) ( Att/j, 2 . 
27rr(l - 2e) 



,2\ e 



(—)5^[2P bb (z,e)G b/P (x 1 /z)G b/P (x 2 ) 



dz 1 — z 
+ (xi <-> sc 2 )l — ( Y <L dx\dx2- 

z z 



(3.11) 



where G bb j P (x) is the bare parton distribution function of 6(6) quark in proton. P bb {z,e) is 
the d-dimensional unregulated {z < 1) splitting function related to the usual Altarelli-Parisi 
splitting kernel ^Hj. P bb (z,e) can be written explicitly as 



Pbb(z,e) 
Pbb(z) 



= Pi 



\ b (z) + ePl b (z), 

n 



= C F 



1 + z 2 



III. .3 Real Corrections from subprocesses 



-C F (l - z). 

-> H+H- + 6(6) 



(3.12) 



In addition to the real gluon emission subprocess bb — > H + H~ + g, there are also subprocesses 
gb(b) — > H + H~~ + 6(6) at this order of perturbation theory, as shown in Fig. 5. 

The contributions from these processes only contain the initial state collinear singularities. 
Using the method described above, we split the phase space into two regions: collinear region 
and non-collinear region. 



a^gbib) -» H+H' + 6(6)) = a^ (gb(b) -» + 6(6)) + a 6 G '( ff 6(6) -> H+fT + 6(6)) 



Ci 



(3.13) 
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Figure 5: The Feynman diagrams for gb(b) —> H + H + b(b). 

Also aP is finite and can be evaluated in four dimensions using standard Monte Carlo tech- 
niques. The differential cross section da^ can be written as 



da 



C 



da 1 



a s r(l - e) / 47T//; 



27rr(l - 2e) V « 



,2\ e 



(--)5- e [2P b ,(z,e)G s/P (x 1 /z)G 6/P (x 2 ) 



(3.14) 



with 



P bg (z,e) = P bg {z) + ePi g {z) } 
P b9 {*) = \[z 2 + (l-z)\ Pi g (z) 



-z(l-z). 



(3.15) 



III..4 NLO QCD Corrected Cross Section for pp -»• bb -> + X 



After adding the renormalized virtual corrections and the real corrections, the partonic cross 
sections still contain the collinear divergences which can be absorbed into the redefinition of 
the distribution functions at NLO. Using the MS scheme, the scale dependent NLO parton 
distribution functions are given as ^ 



Gi/p(x,n f ) = Gi/ P (x) + 



a s T(l -e) Airtf 



2vr T(l - 2e) I fij 



1 dz 



P ij (z)G j/P (x/z). 



(3.16) 



By using above definition, we get a NLO QCD parton distribution function counter-terms 
which are combined with the collinear contributions (Ea. (|3.11j) and Eq. (|3.14l0 to result in the 
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0(a s ) expression for the remaining collinear contributions: 



da 



coll 



da 



a s T(l — e) / 47T/x^ 
27rr(l - 2e) \~~T 



2\ e l 



{2G b /p(xi : fif)Gi /P (x 2 ,fif) 



+ (aci <-> £ 2 )}<£ci<£c2, 



(3.17) 



where 



and 



with 



Af{b^bg) = C F (2ln5 s + 3/2), ^ c = Afln(— ), 



G b / P (x,fj,f) = / —G j/P (x/y,Hf)P bj (y), 
j=b,g Jx y 



P {j (y) = P id ln(5 c i-^-^) - 



(3.18) 



(3.19) 



(3.20) 



y \i s 

We can observe that the sum of the soft (Eq. (|3.9|) ). collinear (Eq . (j3, 1 7 j) ) . and ultraviolet 
renormalized virtual correction (Eq. (|3.3|) ) terms is finite, i.e., 



A% + A$=0, 

Af + A\ + 2A{ c {b -► 65) = 0. 



(3.21) 



The final result for the total O(o s ) correction can be written as the sum of two terms: a 
two-body term a^ and a three-body term a^ . 



a 



(2) 



/ dxidx2da°{G b / P (xi, fif)Gi/ P {x 2 ,i* f )[Ai + A% + 2A s c (b^bg)] 
+ 2G b /p(x 1 , H if)G b /p(x 2 ,Vf) + (xi x 2 )}. (3.22) 



And 



o- (3) = y dxidx 2 [G b/P (x 1 , /j,f)G b/P (x 2 , fif) + (a?i <-> x 2 )]d<J (3) , (3.23) 
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with 



!_Y\Mz{bb -» H+H-g)\ 2 dT 3 + 2 / V|M 3 ( ff 6 -> i^+i^-fo)! 2 ^}. 




(3.24) 



Finally, the NLO total cross section for — > 66 — > H + H + X is given as 



<7 



NLO 



(3) 



(3.25) 



IV. Numerical Results and Discussion 

In this section, we present the numerical results of the cross section for the charged Higgs 
boson pair production via bottom quark fusion at the LHC. In the numerical evaluation, 
we take the SM parameters as: rat = 178.1 GeV, mz = 91.1876 GeV, my/ = 80.425 GeV 
and OiEwim-w) = 1/128. We use the two-loop evolution of the strong coupling a s (fj, r ) with 
a s (triz) = 0.1187 \T3\. We use the CTEQ6L1 parton distribution function for the LO cross 
sections and CTEQ6M for NLO results [20]. The factorization scale is taken as fit = m#±/4 
and the renormalization scale is set to be \i r = mu± by default unless otherwise stated. We 
present the results involving the 'pure' SM-like QCD and the total SUSY QCD corrections 
in following two subsections separately. 

IV.. 1 Results Including Only the SM-like QCD Corrections 

In this subsection, we present the cross sections including only the SM-like QCD corrections. 
It means we consider this process in an Two-Higgs-Doublet model without taking SUSY 
particles into account. When we calculate the virtual one- loop corrections, we only include 
the SM-like QCD one-loop diagrams shown in Fig. 2 and set the renormalization constants 
in Ea. (J3.1|) to be their QCD parts. The MS bottom quark mass mb(^ r ) can be evaluated 
by using the one-loop or two-loop renormalization group improved formula with the bottom 
quark pole mass taken to be mj = 4.7 GeV. They are expressed as 



3 7r |_a s (m&) 



I (co/M 



{m b (n r )) 



QCD 

11 




(4.1) 



12 



(m b (fi r ))% CD 



where 



Aa s (m h ) 
"W ~ o — ~ — ) 

O 7T 



_a s (m b ) 



(co/6o) 



1 + t^(ci - 61) (a s (fj, r 
00 



1 z 11 „ 2 , 1 

47T 3 3 J 7T 

1 51iV- 19n f 1 , 

b i = 7; — 7TT7 — jtA ci = — (1017V- 10n/). 
2vr UN -2n f ' 72vr v /; 



a a (m b )) 

(4.2) 

(4.3) 
(4.4) 



where iV(= 3) is the number of colors and n/(= 5) is the number of active light flavors. In 
our calculation we use (m b (/x r ))^ C '' D to evaluate the LO cross sections and (m(,(/z r ))^ D for 
the NLO cross sections. 
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Figure 6: The dependence of the cross sections on the cutoff 5 S 

Fig.6 shows that our NLO QCD result does not depend on the arbitrary cutoffs S s and 
5 C of the two cutoff phase space slicing method. The two-body(a^ 2 ^) and three-body (cr^) 
contributions and the NLO cross section (cr NLO ) are shown as a function of the soft cutoff 6 S 
with the collinear cutoff S c = <5 s /50. tan/3 = 40, m H ± = 180 GeV and jif = fi r = m H ±. We 
can see the NLO cross section 

a NLO ig 

independent of the cutoffs. In the following numerical 
calculations, we take 6 S = 10" 4 and S c = S s /50. 
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Figure 7: The dependence of the cross sections on factorization scale /x/ and renormalization 
scale fj, r . 

In Fig.7, we show the dependence of the total cross section for pp — ► bb — ► H + H~ + X on 
the renormalization scale and the factorization scale with tan f3 = 40, m H ± = 180 GeV. In the 
left plot of Fig.7, the renormalization scale is taken to be fi r = m,jj± while the factorization 
scale varies in the region Q.lm H ± ~ 3m H ±. In the right plot of Fig.7, we take fif = m H ± and 
0.lm H ± < /i r < 3m H ±. A significant reduction of the scale dependence for the NLO cross 
sections can be observed, thus the reliability of the NLO QCD predictions has been improved 
substantially. 

Fig.8 shows the dependence of the LO and NLO cross sections for pp — > bb — > H + H~ + X 
on the charged Higgs mass m H ±. The values of tan/3 are taken to be 40 (upper lines) 
and 20 (lower lines). m H ± varies from 180 GeV to 500 GeV. The total cross section can 
reach 10 fb for small values of m H ± with very large tan/3. In Fig.9, the dependence of the 
pp — ► bb — > H + H~ +X cross sections on tan/3 are studied for m H ± = 180 GeV and 400 GeV. 
The cross sections increase rapidly with the increment of tan (5 from 10 to 50. 
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Figure 8: The dependence of the cross sections on m H ± . 




Figure 9: The dependence of the cross sections on tan/3. 
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IV.. 2 Results Including the Total SUSY QCD Corrections 

In this subsection, we present the cross sections including all the NLO SUSY QCD corrections. 
The relevant SUSY parameters in our calculation are: the parameters Mq q- ^ and A^ in 
squark mass matrices, the higgsino mass parameter [i and the mass of the gluino m§. The 
squark mass matrix is defined as 

M\ = ( < ^ q ) (4.5) 

9 V a <l m Q m q R J 

with 



mf,, — M? + m 2 , + m 2 z cos 2(3 (1% — e q sin 2 Oyy), 



m 5 - H = ^ + m q + m z cos 2(3e q sin Q w 

a q = A q — //{cot (3, tan/3}, (4-6) 

for {up, down} type squarks. I| and e q are the third component of the weak isospin and the 
electric charge of the quark q. The chiral states qi and qn are transformed into the mass 
eigenstates q\ and q^. 



92 / V 9fi 

Then the mass eigenvalues and m q2 are given by 

("J 1 m °| ) = R^M 2 q {R^ (4.8) 

For simplicity, we assume Mq = = = A t = A b = m g = Msusy = 500 GeV, fi = 200 
GeV. 

In the MSSM, the counter term of can be very large due to the 'pure' SUSY QCD 
(gluino-mediated) diagram for large values of tan (3. The gluino-mediated contributions can 
be absorbed into the tree-level Yukawa couplings )21| . In such a way we obtain the MS bottom 
quark mass mj including the total SUSY QCD contributions, 

(rn b (n, r )fQ CD = (M b (u, r )f CD - ^C F m b {j2[B{ m - ^ sin(2^)(-l)^ m ]K, mj, m? )} 

i=i 

(4.9) 
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where the notations B( tn and Bq IU denote the operations of taking the finite parts of the 
two-point integral functions. The \i r dependence of B functions cancels after summing over 
the sbottom index i = 1,2. As shown in a series papers, the SUSY QCD contributions to the 
bottom quark running mass can be written as |22| 



(m h (u )) SQCD ~ (mb( ^ )} 
(m 6 (M) - 1 + Amb 



QCD 



(4.10) 



where 



2a 



with 



1 a 2 6 2 c 2 

I(a, 6, C) = — -2-T72 2^ lQ g 12 + 1)2 ° 2 l0 § ~~2 + ^ 1o S ( 412 ) 

(a z — o z )(o z — c z )(a z — <r) o z c z a z 

Indeed, if all supersymmetry breaking mass parameters and /x are of equal size, one get 
an interesting limit of Ea. (|4.1H) 23 , 

A • / ^,a s (Q = Msusy) , a [aiq\ 
Ami, = sign{fi) tanp. (4.13) 

37T 

We can see Ami, does not decouple in the limit of large values of the supersymmetry breaking 
masses. The sign of fi is the decisive factor in determining whether the 'pure' SUSY QCD 
corrections will enhance or suppress the cross section for the process of pp — > bb — > H + H~ +X . 

In our calculations we use Ea. ()4.9|) to calculate the MS running bottom quark mass. The 
large SUSY QCD corrections are absorbed in the tree-level bottom Yukawa coupling, and we 
use the bottom quark mass counter term defined in Ea. ()3.1|) to avoid double counting. In 
Fig. 10, we show the dependence of the LO and NLO cross sections on the charged Higgs mass 
mu± which is same as Fig. 8 but including the total SUSY QCD corrections. The curves in 
Fig. 10 are much lower than those in Fig. 8, because after including the 'pure' SUSY QCD 
contributions (mj()v)) is much smaller than (rub^r)) ■, assuming the sign of is 

positive and tan/3 is large (see Ea. (j4.13|0 . In Fig. 11, we plot the tan/3 dependence of the 
cross section including the total SUSY QCD contributions. 
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Figure 10: The dependence of the cross sections on m H ± including the total SUSY QCD 
corrections. 
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Figure 11: The dependence of the cross sections on tan/3 including the total SUSY QCD 
corrections. 
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In summary, we have studied the production of the charged Higgs boson pair via the 
bottom quark fusion in the MSSM including the NLO QCD contributions at the LHC. With 
very large values of tan f3, bb — > H + H~ subprocess can become the dominant mechanism 
in the charged Higgs boson pair production at the LHC. The numerical results of the cross 
sections show that the NLO SUSY QCD corrections are generally significant. We find also 
that the MS bottom quark mass will receive large corrections from the 'pure' SUSY QCD 
contributions which will significantly suppress or enhance the cross section depending on the 
sign of the higgsino mass parameter jjl. The NLO QCD corrections can significantly reduce 
the dependence of the cross sections on the renormalization and factorization scales. 
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